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Abstract 



From Dalitz-plot analyses of 7r^7r+7r+ and £>+ K^n^n^ 

decays, we find evidence for light and broad scalar resonances cr(500) 
^ I and k(800). From a Dalitz-plot analysis of — » 7r~7r+7r+ decays, we 

measure the masses and decay widths of the scalar resonances /q(980) 



Dh! and /o(1370). 



1 Introduction 



The constituent quark model of QCD provides a successful descrip- 
■ tion of pseudo-scalar and vector mesons, which can be conveniently 

arranged into SU (3)^ nonet representations. One also expects a corre- 
sponding scalar meson nonet to exist, and much experimental effort has 
been directed towards identifying its members. The following states 
have been identified thus far, listed by isospin: 

. / = 0: /o(600) or a(500), /o(980), /o(1370), /o(1500), /o(1710) 

• 7=1/2: k(800), K^{U30) 

• 7=1: ao(980), ao(1450) 



*The E791 Collaboration is: CBPF (Brazil), UC Santa Cruz, Cincinnati, CIN VESTAS 
(Mexico), Fermilab, IIT (Chicago), Kansas State, Massachusetts, Mississippi, Princeton, 
Puebla (Mexico), South Carolina, Stanford, Tel Aviv, Tufts, Wisconsin, and Yale. 
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Scalar mesons have traditionally been studied via scattering experi- 
ments. However, in these experiments the mesons can be difficult to 
disentangle from nonrcsonant background due to their broad widths 
and lack of a distinctive angular distribution. Thus some scalar states 
are controversial ((t(500) and k(800)), while others, although accepted, 
have poorly-dctcrmincd parameters (/o(980), aQ(980), and /o(1370)). 
For a recent review, see Ref. [1]. 

Recently, large samples of £>^^ — > 7r~7r+7r+ and — > /r~7r+7r+ 
decays have been used to study scalar mesons (S) via the quasi-two- 
body decay £)+ — > Stt'^, S ^ h~n'^ (charge-conjugate modes are as- 
sumed unless noted otherwise). The scalar masses and decay widths 
are determined by fitting a Dalitz-plot distribution to the square of 
coherently-summed decay amplitudes | X]„^(-D— ^S^Tr)!^. In this pa- 
per, we present such results from Fermilab E791. This experiment 
produced D mesons using a 500 GcV/c tt^ beam incident on Pt and C 
targets. The subsequent D decays were reconstructed using a silicon- 
strip vertex detector, a large-aperture spectrometer, and various par- 
ticle identification detectors. Details of the experiment can be found 
in Ref. [2]. 

2 Dalitz-Plot Formalism 

The analyses select 1?^^ — » 7r~7r+7r+ and — > K~'K~^Tr'^ decays by 
requiring that there be a 3-track vertex well-separated from the pri- 
mary interaction vertex. Pions and kaons are identified using infor- 
mation from two threshold Cerenkov counters. The £>+ — » if ~7r+7r+ 
sample contains 15090 events with about 6% background, while the 
£)+ — »• 7r~7r+7r"'" and Df 7r~7r+7r+ samples contain 1686 and 937 
events, respectively, with about 30% background. The Ktttt sample is 
significantly larger and has a lower level of background because this 
decay is Cabibbo-favored. 

To study the resonance structure of each sample, an unbinned 
maximum likelihood fit is performed. The likelihood function is £ = 
rie-uents + '^sl' where Vg is the probabiHty density function (PDF) 
for signal and is the PDF for background. The latter PDF is ob- 
tained from a fit to events in mass sidebands or from Monte Carlo 
simulation. The signal PDF is Vg = (l/A^^) ^(M) e(m?2, mfs) |^|^ 
where Ng is a normalization factor, g{M) describes the signal shape 
in the tttttt or Kittt mass spectrum, and e{m\2i "^13) is the acceptance 
over the Dalitz plot, including smearing. The amplitude A is the co- 
herent sum of a uniform non-resonant (NR) amplitude and amplitudes 
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of resonances: 



N 

A = aoe**o A + J2 a^e'*" A(mf2, mfs) . (1) 

The coefficients a„ and relative phases (5„ are determined from the 
fit. The amphtude A„ = F^^^ F^^ Mi^^ BW„, where F^^^ and F^^ 
are Blatt-Weisskopf factors that depend on the spin (J) and radii (r^^ 
and Tfi) of the parent D and intermediate resonance; Air/"* is tlie 
angular-momentum-conserving factor (—2)'^ p"^ q"^ Pj {cos 9), where p, q 
are the momenta of the D and one of the daughters of the intermediate 
resonance in the resonance rest frame, and cos^ = P ■ Q', and BW^ 
is the relativistic Breit-Wigner propagator [m^ — — im„r„(?7?,)] , 
where m is the invariant mass of the track pair forming a resonance 
{mi2 or rrii^), rn„ is the resonance mass, and r„(m) is the mass- 
dependent width. More details of these expressions can be found in 
Ref. [3]. Finally, each resonant amplitude is Bose-symmetrized with 
respect to the identical pions: A„ = ^[(12)3] + ^[(13)2]. 

3 Results 

The results of the likelihood fits are given in Tables 1-3. The tables list 

the coefficients a„ and phases S,^ rcsiilting from the fits, along with the 
decay fraction for each amplitude. The decay fraction for an amplitude 
is defined as its intensity integrated over the Dalitz plot, divided by the 
integrated intensity of all amplitudes coherently summed. To assess the 
quality of the fit, a fast Monte Carlo program was developed that pro- 
duces binned Dalitz-plot densities according to signal and background 
PDFs, including detector smearing effects. The difference between this 
density distribution and that of the data - summed over all bins - is 
taken as the of the fit. The number of bins is the number of degrees 
of freedom. 

Table 1 gives the results for D+ — »7r~7r"'"7r"'". The left- most column 
lists results obtained when fitting to established resonances, while the 
right-most column lists results obtained when an additional scalar res- 
onance is included. For the latter fit, the per degree of freedom (^) 
is significantly improved, and the non-resonant fraction constitutes a 
much smaller component of the total width (in the left-most column 
it dominates). The mass and width obtained by the fit for the ad- 
ditional resonance (denoted "cr") are m = 478^23 17 MeV/c^ and 
F = 324^40 ± 21 MeV/c^. As a test wc also fit the data with vector, 
tensor, and toy (Breit-Wigner with no phase variation) models for the 
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additional amplitude. In all cases we obtain a higher ^ jv than that 
obtained with the additional scalar amplitude. 

Tabic 2 gives the results for ^7r^7r+7r+. The left-most column 
lists results obtained when fitting to established resonances, while the 
other two columns list results obtained when the /9°(770) and /9°(1450) 
are cxchidcd. The parameters determined from the fit are the mass 
and width of the /"(ISTO), and the mass and coefficients g^^ 5^ of 
a coupled-channel Breit-Wigner width [4] for the /°(980). The results 
aremj. (^g^o) = 1434±18±9 MeV/c^, (^g^^^ = 172±32±6 MeV/c^, 

and ^g^^^ = 977 ± 3 ± 2 MeV/c^, = 0.09 ± 0.01 ± 0.01, gj^ = 
0.02 ± 0.04 ± 0.03. We also fit the data using the same Breit-Wigner 
for the /q(980) as that used for the other resonances; the resulting fit 
is nearly as good and the decay fractions and phases are essentially 
unchanged. The isoscalar plus tt^ components correspond to over 90% 
of the _D+ ^7r~7r"'"7r"'" decay width. 

Tabic 3 gives the results for ti^ . The left-most column 

lists results when fitting only to established resonances and fixing the 
parameters of the scalar /Cq (1430). The middle column lists results ob- 
tained when the ifQ(1430) parameters are allowed to float, and, in ad- 
dition, Gaussian form factors Fj^-F^ = exp(— p^r|j/12)-exp(— p^r|j/12) 
are used for the Kq{1430) to account for the size of the decaying 
mesons. The meson radii and (also appearing in Blatt-Weisskopf 
factors) are now determined from the fit. The right-most column lists 
results when an additional scalar resonance is included (also with Gaus- 
sian form factors in its amplitude). This last fit yields a significantly 
improved x^/'^j and the NR fraction is reduced from about 90% to 13%, 
which is more consistent with expectations. The mass and width of the 
additional resonance (denoted "/t") are m = 797 ±19 ±43 MeV/c^ and 
r = 410 ± 43 ± 87 MeV/c^. Including this additional resonance, the fit 
obtains Kq{1430) parameters Jri^.^^^gg^ = 1459 ± 7 ± 5 MeV/c^ and 
^K*{iA3Q) = 175±12±12 MeV/c^: these values are shifted by +47 MeV 
and —119 MeV, respectively, with respect to PDG values [6]. The me- 
son radii are =5.0 ± 0.5 GeV"^ and r^ = 1.6 ±1.3 GeV'^ We also 
fit the data with vector, tensor, and toy (Breit-Wigner with no phase 
variation) models for the additional amplitude; in all cases we obtain 
a higher x^/^^ than that obtained with the k amplitude. 

In summary, we have performed Dalitz-plot analyses for the de- 
cays D+ — » 7r~7r"'"7r+, Df 7r~7r"'"7r+, and D+ — > if ~7r+7r+. For the 
first sample, including an additional scalar resonance with mass near 
500 MeV/c^ significantly improves the quality of the fit. For the £>+ 
K~Tr~^7r^ sample, including a scalar resonance with m « 800 MeV/c^ 
improves the quality of the fit. When this resonance is included, the fit 
obtains a mass and decay width for the i^Q (1430) that are shifted with 
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Estab. Resonances 


Adding cr(500) 




Fraction (%) 


Fraction (%) 




Phase 


Phase 


NR 


38.6 ±9.7 


7.8 ±6.0 ±2.7 




(150 ± 12)° 


(57.3 ±19.5 ±5.7)° 


fTTT"'" 


— 


46.3 ±9.0 ±2.1 






(206 ±8.0 ±5.2)° 


p'J(770)7r+ 


20.8 ±2.4 


33.6 ±3.2 ±2.2 




0° (fixed) 


0° (fixed) 


/o(980)7r+ 


7.4 ± 1.4 


6.2 ± 1.3 ±0.4 




(152 ±16)° 


(165 ±11 ±3.4)° 


/2(1270)^+ 


6.3 ±1.9 


19.4 ±2.5 ±0.4 




(103 ± 16)° 


(57.3 ±7.5 ±2.9)° 


/o(1370)7r+ 


10.7 ±3.1 


2.3 ± 1.5 ±0.8 




(143 ±9.7)° 


(105 ±18 ±0.6)° 


^•^(1450)77+ 


22.6 ±3.7 


0.7±0.7±0.3 




(46± 15)° 


(319 ± 39 ± 11)° 




254/162 


138/162 



Table 1: Results of the maximum likelihood fit to the -^ir n'^n'^ Dalitz 
plot (from Ref. [3]). 

respect to PDG values. For the Z)+ ^7r~7r+7r+ sample, the Dalitz-plot 
fit yields new values for the masses and decay widths of the /q(980) 
and /o(1370) scalar mesons. 
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Estab, Resonances 


Wo p {^i iyj) 


INO p ^14oU ) 




Fraction (%) 


Fraction (%) 


Fraction (%) 




Magnitude 


Magnitude 


Magnitude 




Phase 


Phase 


Phase 


NR 


0.5 ±1.4 ±1.7 


7.5 ±4.8 


5.0 ±3.8 




0.09 ±0.14 ±0.04 


0.36 ±0.12 


0.30 ±0.12 




(181 ±94 ±51)° 


(165 ±23)° 


(149 ± 25)° 


/O (980)77+ 


56.5 ±4.3 ±4.7 


58.0 ±4.9 


54.1 ±4.0 


1 (fixed) 


1 (fixed) 


1 (fixed) 




0° (fixed) 


0° (fixed) 


0° (fixed) 


p"(770)7r+ 


5.8 ±2.3 ±3.7 




11.1 ±2.5 




0.32 ±0.07 ±0.19 




0.45 ± 0.06 




(109 ±24 ±5)° 




(81 ±15)° 


/2(1270)7r+ 


19.7±3.3±0.6 


22.2 ±3.3 


20.8 ±3.0 




0.59 ±0.06 ±0.02 


0.62 ±0.06 


0.62 ±0.05 




(133 ± 13 ±28)° 


(109 ± 11)° 


(124 ± 11)° 


/o(1370)7r+ 


32.4 ±7.7 ±1.9 


30.4 ±6.9 


34.7 ±7.2 




0.76 ±0.11 ±0.03 


0.72 ±0.11 


0.80 ±0.11 




(198 ± 19 ± 27)° 


(1.5(j ± 19)" 


(159 ± 14)" 


p^(1450)7r+ 


4.4 ±2.1 ±0.2 


5.8 ± 2.2 






0.28 ±0.07 ±0.01 


0.32 ± 0.06 






(162 ± 2(j ± 17)" 


(144± 20)" 






72/68 


93/68 


104/68 



Table 2: Results of the maximum likelihood fit to the Df ^tt tt'^tt'^ Dalitz 
plot (from Ref. [4]). 
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Estab. Resonances 


Floating K*(1430) with 


Adding k(800) 






Gaussian form factors 






Fraction (%) 


Fraction (%) 


Fraction (%) 




Magnitude 


Magnitude 


Magnitude 




Phase 


Phase 


Phase 


NR 


90.9 ±2.6 


89.5 ± 16.1 


13.0 ±5.8 ±4.4 




1.0 (fixed) 


2.72 ±0.55 


1.03 ±0.30 ±0.16 




0° (fixed) 


(-49 ±3)° 


(-11 ±14 ±8)° 








47.8 ±12.1 ±5.3 








1.97± 0.35 ±0.11 








(187±8± 18)° 


/r*(892)7r+ 


13.8 ±0.5 


12.1 ±3.3 


12.3 ± 1.0 ±0.9 




0.39 ±0.01 


1.0 (fixed) 


1.0 (fixed) 




(54 ±2)° 


0° (fixed) 


0° (fixed) 


Jfo* (1430)77+ 


30.6 ±1.6 


28.7 ± 10.2 


12.5 ± 1.4 ±0.5 




0.58 ±0.01 


1.54 ±0.75 


1.01 ±0.10 ±0.08 




(54 ±2)° 


(6± 12)° 


(48±7± 10)° 


Jfl (1430)77+ 


0.4 ±0.1 


0.5 ±0.3 


0.5 ±0.1 ±0.2 




0.07 ±0.01 


0.21 ±0.18 


0.20 ±0.05 ±0.04 




(33 ± 8)° 


(-3 ±26)° 


(-54±8±7)° 


iif* (1680)77+ 


3.2 ±0.3 


3.7 ± 1.9 


2.5 ±0.7 ±0.3 




0.19 ±0.01 


0.56 ±0.48 


0.45 ±0.16 ±0.02 




(66 ± 3)° 


(36 ± 25)° 


(28 ± 13 ± 15)° 




167/63 


126/63 


46/63 



Table 3: Results of the maximum likelihood fit to the 7r"'"7r^ Dalitz 

plot (from Ref. [5]). 
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